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Abstract 

Let G be a connected planar (but not yet embedded) graph and F a set of additional edges 
not yet in G. The multiple edge insertion problem (MEI) asks for a drawing oi G + F with 
the minimum number of pairwise edge crossings, such that the subdrawing of G is plane. An 
optimal solution to this problem approximates the crossing number of the graph G + F. 

Finding an exact solution to MEI is NP-hard for general F, but linear time solvable for the 
special case of |F| = 1 (SODA 01, Algorithmica) or when all of F are incident to a new vertex 
(SODA 09). 

The complexity for general F but with constant k = \F\ was open, but algorithms both with 
relative and absolute approximation guarantees have been presented (SODA 11, ICALP 11). 
We show that the problem is fixed parameter tractable (FPT) in k for biconnected G, or if the 
cut vertices of G have degrees bounded by a constant. We give the first exact algorithm for this 
problem; it requires only 0{\V{G)\) time for any constant k. 


*P. Hlineny has been supported by the Czech Science Foundation project 14-03501S. 
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1 Introduction 


The crossing number cr(G) of a graph G is the minimum number of pairwise edge crossings in a 
drawing of G in the plane. Finding the crossing number of a graph is one of the most prominent 
combinatorial optimization problems in graph theory and is NP-hard already in very restricted 
cases, e.g., even when considering a planar graph with one added edge (cf. the MEI problem 
for k = 1 later). The problem has been vividly investigated for over 60 years, but there is still 
surprisingly little known about it; see e.g. 30 for an extensive reference. While in general, there 


exists a c > 1 such that the crossing number cannot be approximated within a factor c in polynomial 
time [^, several approximation algorithms arose for special graph classes. 

For general graphs with bounded degree, there is an algorithm that approximates the quantity 



In this paper, we follow the latter idea and concentrate on the Multiple Edge Insertion problem 
MEI(G, F), to be formally defined in the section hereafter. Intuitively, we are given a planar graph 
G, and ask for the best way (in terms of total crossing number) to planarly draw G and insert a 
set of new edges F into G such that the final drawing of G + F (i.e., of the graph including the 
new edges of F) restricted to G remains planar. 

This problem is polynomial-time solvable for |E| = 1 and in the case when all edges of F 
are incident to a common vertex [^, but NP-hard for general F 33 . Moreover, an exact or 


at least approximate MEI solution constitutes an approximation for the crossing number of the 
graph G + F [^. Considering constant k := \F\, there have been two different approximation 
approaches |13| and [10] ; the former one directly targets the crossing number and achieves only a 
relative approximation guarantee for MEI; the latter one first specifically attains an approximation 


of MEI with only an additive error term, and then uses 11 to deduce a crossing number approx- 

in fact turns out 


imation. While the former one is not directly practical, the latter algorithm [^ 
to be one of the best choices to obtain strong upper bounds in practice [^. 

In this paper, we show that for every constant k and under mild connectivity assumptions, there 
is an exact linear time algorithm, which has so far been an open problem even for k = 2. In terms 
of parameterized complexity, our algorithm is in EFT with the parameter k = |E|. 

Theorem 1. Let G be a planar connected graph on n vertices, and F a set of k new edges (vertex 
pairs, in fact) where k is a constant. If G is biconnected, or the maximum degree of the cut vertices 
of G is bounded by a constant, then the problem MEI(G, E) is solvable in 0{n) time. 


We also mention that while the crossing number itself is in FPT w.r.t. the objective value 
already a planar graph with one added edge may have unbounded crossing number. 


ISpG 


Both the aforementioned absolute MEI-approximation [1C)| and our new approach can use 11 to 


obtain the same relative ratio for approximating the crossing number. However, our new approach 
does so without any additional additive term: 


Corollary 2. Using the theorem relating an optimum MEI(G, E) solution to the crossing number 
of the graph G + F JIT] /, Theorem 0 gives a polynomial time kA-approximation for the crossing 
number of G + F with constant k = |E|, where G is a planar graph and A is its maximum degree. 
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Organization. After formally defining our setting in the next section, we will concentrate on 
the still NP-hard problem Rigid MEI in Section]^ i.e., MEI under the restriction that the planar 
embedding of G is fixed. In Section]^ this algorithm is at the core of a dynamic programming over 
a decomposition tree of G, in order to obtain an FPT algorithm for the general MEI, i.e., when 
any planar embedding of G is allowed. The latter constitutes the result for Theorem 

2 Preliminaries 

We use the standard terminology of graph theory. By default, we use the term graph to refer to a 
loopless multigraph, i.e., we allow parallel edges but no self-loops. If there is no danger of confusion, 
we denote an edge with the ends u and v chiefly by uv. 

A drawing of a graph G = {V, E) is a mapping of the vertices V to distinct points on a surface 
S, and of the edges E to simple (polygonal) curves on S, connecting their respective end points 
but not containing any other vertex point. Unless explicitly specified, we will always assume S 
to be the plane (or, equivalently, the sphere). A crossing is a common point of two distinct edge 
curves, other than their common end point. Then, a drawing is plane if there are no crossings. 
Plane embeddings form equivalence classes over plane drawings, in that they only define the cyclic 
order of the edges around their incident vertices (and, if desired, the choice of the outer, infinite 
face). A planar graph is one that allows a plane embedding. A plane graph is an embedded graph, 
i.e., a planar graph together with a planar embedding. 

Given a drawing D of G, let ct{D) denote the number of pairwise edge crossings in D. The 
crossing number problem asks for a drawing D° of a given graph G with the least possible number 
ct{D°) =: cr(G). By saying “pairwise edge crossings” we emphasize that we count a crossing point 
X separately for every pair of edges meeting in x (e.g., i edges meeting in x give ( 2 ) crossings). 

Definition 3 (Multiple edge insertion, MEI and rigid MEI). 

Consider a planar, connected graph G and a set of edges (vertex pairs, in fact) F not in E{G). We 
denote hy G + F the graph obtained by adding F to the edge set of G. 

Let Go be a planar embedding of G. The rigid multiple edge insertion problem r-MEI(Go,T) 
is to find a drawing D of the graph G + F with minimal ct{D) such that the restriction of D to G 
is the plane embedding Gq. The attained number of crossings is denoted by r-ins(Go, F). 

The multiple edge insertion problem MEI(G, F) is to find an embedding Gi of G (together 
with the subsequent drawing D as above), for which r-MEI(Gi,T) attains the minimum number 
of crossings. The latter is denoted by ins(G, F). o 

Herein, we will also deal with the weighted crossing number, i.e., we have edge weights w : E{G) —)• 
N+ U {00}, and a crossing between two edges 61,62 accounts for the amount of 66(61) • w^e-z) in 
the above crossing functions. Specially, for the MEI problem variants, we shall consider integer 
weights on the edges of G but not on F (i.e., the weight on F is always 1). Although this is not a 
noteworthy strengthening of Theoremby itself, the weights on E{G) will be useful in the recursive 
processing of the non-rigid case, cf. Section 

Given a plane embedding Gq of G, we define its (geometric) dual Gg as the embedded multigraph 
that has a (dual) vertex for each face in Gq; dual vertices are joined by a (dual) edge for each (primal) 
edge shared by their respective (primal) faces. The weight of a primal edge gives rise to the length 
(of same value) of its dual edge. The cyclic order of the (dual) edges around any common incident 
(dual) vertex v*, is induced by the cyclic order of the (primal) edges around the (primal) face 
corresponding to v*. 
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We refer to a path/walk in Gq as to a dual path/walk in Go, and we speak about a dual path/walk 
TT in Go between vertices u, v if the tt starts in a face incident with u and ends in a face incident 
with V. We shortly say a route from u to v {a u-v route) to mean a dual walk between vertices u, v. 

For any drawing D, let cr£)(X, y) denote the number of crossings between edges of X and 
edges of y in D, and let cid{X) := ciniX, X). It is well established that the search for an optimal 
solution to the crossing number problem can be restricted to so-called good drawings: any pair of 
edges crosses at most once, adjacent edges do not cross, and there is no point that is a crossing of 
three or more edges. A simple extension of this finding to the setting of MEI is presented next, in 
Lemma 01 

The following technical results will be used to restrict how “complicated” drawings of the edges 
of F may look in an optimal solution of a MEI(G, F) or r-MEI(G, F) instance. Note that, although 
both the claims are formulated for the rigid version, they easily imply the same for the ordinary 
(non-rigid) MEI problem. 

Lemma 4. Consider a (weighted) instance r-MEI(G, E) of the rigid MEI problem. In any optimal 
solution o/r-MEI(G, E), any two edges of F eross at most once, and they have no erossing if they 
share a eommon endvertex. Moreover, if the weights of the edges in F equal 1 and there exists a 
drawing D of G + F such that ctd{E{G), F) = c, then r-ins(G, E) < c -|- . 

Proof. The proof simply repeats, for this special case of rigid MEI, the folklore “arc exchange” 
argument from the crossing number theory. Eor the second claim, we observe that since the edge 
weights in E are all 1, it is ctd{F,F) < □ 

Corollary 5. Consider a (weighted) instance r-MEI(G, E) such that the weights of the edges in 
F equal 1, and let f £ F. Assume that Di and D 2 are two drawings 0 / G -|- E sueh that Di — f 
is identical to D 2 — f, and that cr£)^(E(G), {/}) — cr£)^{E{G),{f}) > Then cr(Ei) > 

r-ins(G,E), i.e., Di is not an optimal solution o/r-MEI(G, E). 

This claim might look rather weak at first sight, with respect to the required large difference 
d := cr£)j(/,E(G)) — cr/^^C/,E(G)). However, one can actually easily construct examples in which 
d = H(|E|) and yet Di is an optimal solution to r-MEI(G,E). 

Proof. Let E = E{G), k = |E| and F' = F \ {/}. Using Lemma 0 we estimate 

r-ins(G, E) < cr^i^ (E, E) = cr^^ (E, E') -b cr^^ (E, {/}) -b 

= [ crui (E, E') -b croi (E, {/}) ] - cvd^ (E, {/}) -b cin^ (E, {/}) -b 

< cr£)j(E, E)-b 0 < cr(Ei). 

□ 

3 Rigid MEI 

In this section we give an EPT algorithm for solving the rigid version r-MEI(G, E), parameterized 
by A: = |E|. G is hence a plane graph (i.e., with a fixed embedding) throughout this section. Recall 
that the r-MEI(G, E) problem is NP-hard for unrestricted k. 

We first illustrate the simple cases. Solving r-MEI(G, {uv}), the fixed embedding edge insertion 
problem with A: = 1, is trivial. Augment dual G* with edges of length 0 between the terminals 
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u,v (technically, new vertices in G*) and their respective incident faces (vertices in G*), to suit the 
above definition of a u-v route in G. Realizing a route for uv means to draw uv along it within G. 
If the shortest route has length I, realizing it attains r-ins(G, {urc}) = the smallest number of 
crossings in the rigid MEI setting. 

For k > 2, the situation starts to be more interesting: not every pair of shortest routes gives 
rise to an optimal solution of r-MEI(G, F) since there might arise a crossing between the two edges 
of F. The question, for A: = 2, is whether some pair of shortest routes of the two edges in F can 
avoid crossing each other. Since it is generally not feasible to enumerate all shortest routes, we 
cannot check this by brute force and a more clever approach is needed. Even worse, for larger 
values of k we can encounter situations in which optimal solutions of r-MEI(G, F) draw edges of F 
quite far from their individual shortest routes (in order to avoid crossings with other edges of E). 

On a very high level, our approach to finding a drawing D oi G + F that is an optimal solution 
to r-MEI(G,E), can be described as follows: 

(I) We guess, for each pair /, /' G F, whether / and f will cross each other in D. Since k = |E| 
is a parameter, all the possibilities can be enumerated in FPT time. 

(II) Let X C (^) be a (guessed) set of pairs of edges of F. We find a collection of shortest routes 
for the edges of E in G under the restriction that exactly the pairs in X cross; D is obtained 
by inserting the edges of E along their computed routes. As we will see, we may restrict our 
attention only to routes pairwise crossing at most once. 

(Ill) We select D which minimizes the sum of |X| and of the lengths of the routes found above. 


3.1 Handling path homotopy of rontes 


The core task of the scheme (I) - (III) is to find a collection of shortest routes under the restriction 
that every route avoids crossing certain other routes (note; none of these routes are fixed in advance). 
Although this problem may seem equivalent, in the dual, to the notoriously hard problem of shortest 


disjoint paths in planar graphs 114,28 , this is fortunately not the case since our routes may freely 


share their sections as long as they do not cross. We give a solution of the core task which is greedy 
in the sense that each route of E in G is minimized regardless of the other routes of E. The key to 
this solution is the concept of a path homotopy in the plane. 

In a brief and rather informal topological view, consider the sphere with a finite set of point 
obstacles. Two simple curves a, a' with the same endpoints are homotopic if there exists a homeo- 
morhpism (a continuous deformation) of a to a' that fixes the endpoints and otherwise avoids all 
the obstacles. For example, if a, a' are disjoint except at the common ends, then they are homo¬ 
topic if and only if one of the two open regions bounded by a U a' is obstacle-free. In our case, the 
obstacles are the ends V{F) of the edges of E (as given by the fixed embedding of G), where each 
endpoint is “blown up” into a small open disc. Then, given the homotopy classes hom{a), hom{l3) 
of two curves a,/?, one can decide whether a and /3 are “forced to cross”—although, a and (5 may 
cross if they are not forced to, such an unforced crossing can as well be avoided in our case. 

Instead of the above classical algebraic-topology setting of homotopies, in this paper we prefer 
to deal with path homotopy in a combinatorial setting. This setting is closely inspired by the 
discrete-geometry view of boundary-triangulated 2-manifolds by Hershberger and Snoeyink 21 . 
In the first step, we “triangulate” the point set V{F) (our obstacles) using transversing paths in 
the embedding G. A transversing path between vertices x, y of G is a path whose ends are x, y 
and whose internal vertices subdivide some edges of G. Let T be the union of these transversing 
paths and G' denote the corresponding subdivision of G. In order to avoid a terminology clash with 
graph triangulations, we will call T in the pair {G',T) a trinet of G. Formally (where V{F) = N): 
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Definition 6 (Trinet). Let G be a connected plane graph and N C V{G), |A^| >4. A plane graph 
T such that V{T) n V{G) = N is called a trinet of G if the following holds: 

a) T is a subdivision of a 3-connected plane triangulation on the vertex set N (in particular, 
every face of T is incident with precisely three vertices of iV), and 

b) there exists a subdivision G' of G such that V{G') \ V{G) = V{T) \ N, E{G') n E{T) = 0 
and the union G' U T is a plane embedding. 

The pair {G\T) is a full trinet of G. The vertices in N{T) := N are called trinodes of T, the 
maximal paths in T internally disjoint from N are triedges and their set is denoted by I{T), and 
the faces of T are tricells. Note that the triedges of T are transversing paths of G. We refer to 
Figure for a brief illustration of this definition. o 

Second, we focus on terms related to path homotopy in a full trinet (G', T) of a plane graph G. 
Moreover, while we have implicitly perceived a route of uv in G (i.e., a dual walk from u to v) as 
an arc drawn from u to v, we would also like to describe a topological “alley” for all u-v arcs of a 
similar kind (and same number of crossings) in the embedding G'UT. With it we gain combinatorial 
abstraction and will later be able to avoid unforced crossings with other routes. 

Definition 7 (Alley and T-sequence). Let {G',T) be a full trinet of a plane graph G. Consider a 
route TT between u,v G V (G) in the graph G' U T. Then V (vr) = {cfo, 4>i, ■ ■ ■, 4>m} where each dual 
vertex cfi of vr is an open face of G' U T. Let these faces {cpo, cpi,, (fm) be ordered along vr such 
that 4>o is incident to u and (fm incident to v. Let (ei, 62 ,..., Cm) C E{G' U T) be the sequence 
of the primal edges of the dual edges of vr, ordered from fiQ to 4>m- As a point set, each edge Ci is 
considered without the endpoints. 

a) The union {u, u} U UI^o is called the alley of vr (or, an alley between u, v). 

b) Let (e'^,..., e^) C (ei, 62 ,..., Cm) be the restriction to E(T), and let (pi,p 2 ,... ,Pi) C I(T) be 
the sequence of triedges such that pi contains the edge e( for v = 1,..., ^. Then {pi,p 2 ,... ,P£) 
is called the T-sequence of tt from u to v (or, of the corresponding alley from utov). o 

A route vr crosses a triedge p if the alley of vr contains one of the G'-edges forming p. The 
T-sequence of vr hence describes the unique order (with repetition) in which its alley crosses the 
triedges of T. Usually, we shall consider only the case of tt, u G N{T). 

A route may, in general, cross the same triedge q many times in one place (switching “there and 
back”). However, such a situation may be easily smoothed down to one or no crossing, and this 
can be formalized by the notion of reducing a T-sequence as follows: if 5 = {pi,P 2 , ■ ■ ■ ■,Pe) is a T- 
sequence such that pi = pj+i for some 1 <i < I, then the subsequence S' = {pi,..., Pi-i, Pi +2 ■,■■■, Pi) 
is called a one-step reduction of S. A subsequence S'* C S' is a reduction of S (or S reduces to S* ) 
if S* results from a sequence of one-step reductions of S. 

It comes as no surprise that T-sequences are closely related to the homotopy concept: 

Remark 8 . Consider a trinet T in the sphere. One can show that two arcs with the same fixed 
endpoints are path-homotopic (in the sphere with the obstacles formed by the trinodes of T) if, 
and only if, their T-sequences can be reduced to the same subsequence. However, since we are not 
going to directly use this fact, we refrain from giving this as a formal statement in the short paper. 
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Figure 1: An example of a trinet in a plane graph (see Definition [^ : the underlying plane graph 
G is in black, the trinodes in thick red and the triedges in blue. 



Figure 2: An example (see Definitionj^: the T-sequence of the u-v route depicted in red is {p, q, r, s) 
from u to V. It is a proper T-sequence from u to u (Definition |10[). 
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3.2 T-sequences of potential shortest routes 


Our goal is to look for shortest routes in G of a given homotopy, and we slightly generalize the 
setting to allow for a connected plane graph G with edge weights w: E{G) —)> N+ U {oo}. For a full 
trinet {G\T) of G, we define the edge weights of G' U T as follows: w'{p) := 0 for all p G E{T) 
and w'{e') := w{e) where e' G E{G') is obtained by subdividing e G E{G). This w' is the weight 
induced by w in the trinet {G',T). We give the same weights w' also to the edges of the geometric 
dual of G' U T. If a is the alley of a route vr between vertices x, y in G' U T, then the length of a 
equals the length of vr, i.e., the sum of the tc'-weights of the dual edges of vr. 

With the help of the framework developed in the previous section, we can now give an (again 


informal) high-level refinement of our solution steps (I) - (III) of r-MEI(G, E') as follows; 


(IV) 


(V) 


Consider a trinet T of G on the trinodes V{F). If we fix a (realizable) T-sequence S, then we 
can use established tools, namely an adaptation of the idea of the funnel algorithm j6,29 , to 
efficiently compute a shortest alley among those having the same T-sequence S. For uv € E 
of weight 1, if we compute an alley a between u,v of length d., then we can easily draw the 
new edge uv as an arc in a with d weighted crossings. 

Suppose that, for f = 1, 2, ctj is a shortest alley between Xi and pi having the T-sequence Si. 
Then, as detailed later in Lemma 17 and also Claim 20 we can decide from only 5'i,S'2 


(VI) 


whether there exist arcs from xi to yi in ai and from X 2 to y 2 in a 2 , which do not cross (note 
that n 02 may be nonempty and yet there may exist such a pair of non-crossing arcs). 
Moreover, if the two arcs cross then it should be only once. 

Consequently, it will be enough to loop through all “suitable” T-sequences for every edge 
of F and independently perform the steps (IV)[ 
to get an optimal solution of r-MEI(G, F) as in 


(V) for each combination of them, in order 
The point is to bound the number of 


(III) 


T-sequences that have to be considered, in terms of only the parameter k = |T|. 


We first resolve the last point [(VI)] which is a purely mathematical question. In order to achieve 
the goal, we will build a special trinet of G along (at least locally) shortest dual paths between the 
trinodes in G (Definition]^. Then we will be able to restrict our attention to special T-sequences 
of bounded length (Definition and Lemma [XI] ) . 

Definition 9 (Shortest-spanning trinet). Let (G',T) be a full trinet of a plane graph G, and let 
the weights w' in {G',T) be induced by weights w in G. For a triedge q G /(T), every internal 
vertex t of y is incident with two edges e, e' of G' of weight w'{e) = w'{e') which we call the weight 
of t. The transversing weight of q equals the sum of the weights of the internal vertices of q. 

A triedge q G I{T) between trinodes x,y is locally-shortest if the transversing weight of q is 
equal to the length of a shortest dual path vr in G' U T between x, y, such that vr is contained in(!) 
the union of the two tricells incident with q. Similarly, q is globally-shortest if the transversing 
weight of q is equal to the dual distance between x, y in G' U T. 

We say that T has the shortest-spanning property if every triedge in I{T) is locally-shortest, 
and there exists a subset of triedges J E I (T) forming a connected subgraph of T spanning all the 
trinodes such that every triedge in J is globally-shortest. o 


Definition 10 (Proper T-sequence). Consider a trinet T and trinodes u v € N(T). A nonempty 
sequence S = {pi,P 2 , ■ ■ ■ ,Pm) C /(T) of triedges of T (repetition allowed) is a proper T-sequence 
from u to V the following holds: u is disjoint from pi but there exists a tricell 9q incident with 
both u and pi, u is disjoint from pm but there exists a tricell 9m incident with both v and pm, and 
each two consecutive triedges pi,pi+i are distinct and incident to a common tricell 9i for 1 < z < m. 
Empty 5 is a proper T-sequence from u to v u,v are incident to a common tricell 9q. 


o 




















Recalling that T is a subdivision of a triangulated graph, we immediately get the following; 


Claim 11. For every proper nonempty T-sequence S, the sequence of tricells {Oq^Oi^ ... ,9^ 


Definition 10 is uniquely determined by T and S. 


as in 

n 


For empty proper S, a tricell 6q from Definition 10 is not unique (there are two choices for it). 
However, since any of the two choices of Oq incident with both u, v will work for us in the same way 
in a shortest-spanning trinet, we simply make an arbitrary deterministic choice of 6q for S = 0 and 


extend the scope of Claim also to empty proper T-sequences. 

What follows is the crucial finding that makes our algorithm to work: 


Lemma 12. Consider an instance r-MEI(G, F") where G is a connected plane graph. Let {G',T) 
be a full trinet of G having the shortest-spanning property. There exists a set {vr/ : / G F} where 
'^f for f = uv is a route in G' UT between the trinodes u,v, such that the following hold: 

a) There exists an optimal drawing D of G + F with r-ins(G', F) crossings such that each edge 
f £ F is drawn in the alley of -Kf, and no two edges of F cross each other more than once. 

b) The T-sequence Sf of each nf is a proper T-sequence, and no triedge occurs in Sf more than 
8k^ times where k = |F|. 


Note that optimality of the number r-ins(G, F) ensures that, in a), no two edges of F cross in 
D unless they are forced to (by their given alleys). 


Proof. In the scope of this proof, we shall use the following special terminology and notation. For 
simplicity, we use the symbol / both for an edge / G F and for the arc representing / in a specific 
drawing of G + F (more generally, of (G' + F)UT). We similarly consider a triedge p G I{T) also as 
the arc representing p in G'. If x, y are two points on any arc b, then let b[x, y] denote the section 
of the arc from x to y. 

Consider the given shortest-spanning trinet and the corresponding plane embedded graph G'UT 
with edge weights w' induced by given w of G. We will implicitly assume that every arc a drawn 
in U T avoids crossing the vertices and a intersects G' U T in finitely many points, i.p., the 
embedding and the arcs may be restricted to polygonal lines. For any arc b with the ends u, v, we 
define the T-sequence of b from u to u as the sequence (with repetition) in which b intersects the 
triedges of T. We dehne the transversing weight of b, shortly t-weight, as the sum of the tcCweights 
of the edges of G' U T crossed by b, and denote it by t^fib). 

We choose an optimal drawing F of G^ -|- F which, at the same time, minimizes the combined 
length of the T-sequences of the edges of F, i.e., the number of crossings between F and the 
trinet T. Recall from Lemma any two edges of F cross at most once, and they have no crossing 
if they share a common endvertex. For / G F, let Fj = (pj,... be the T-sequence of / and 

let xj-,... respectively, denote the points at which the arc of / intersects the triedges of T. 

The first task is to prove that each S'/ is a proper T-sequence. 

We start with a stronger technical claim: if, for some / G F and j > i, it is p/ = pj- and the 

simple loop a := p).[x)-,Xy] U /[x/,Xj] is contractible (i.e., with no trinode inside), then we get a 
contradiction to the choice of D above. Indeed, we may assume that / and j > i are chosen such 
that a encloses minimal area in the drawing D. By the minimality of a, no triedge crosses the 
interior of /[xpxj,] twice (all P/^^, • ■ • )P/~^ distinct). However, since the interior enclosed by 
a contains no trinode, the previous implies that no triedge other than p). may intersect a, and so 
j = i + 1. Consequently, since the triedge p). is locally shortest in T, the t-weights of the considered 
section satisfy (p/[xp x^]) < tio'(/[x/, xj-]). If we re-route / closely along p).[x/,xj.] (without 
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crossing p^), then this change does not increase the crossing number by the inequality of t-weights, 
but the T-sequence of / gets shorter (see in Figure]^. Hence, it contradicts our choice of D. 

Now we get back to Sf being a proper T-sequence. If Sf is empty, then the statement is 
trivial. If 5"/ contains consecutive repeated triedge pj = for some 1 < i < mf, then the 
above contradiction directly applies. Assume now that f = uv and the triedge pj. is incident 
with the starting trinode u. Then we can apply the same contradiction to the contractible loop 
a := Xj] U /[u, xj-]. The remaining properties of proper T-sequences follow trivially. 


The last and most difficult step is to prove that no triedge repeats in 5/ too many times, for 
each f G F. Again, if p^j = p’j = p' for i / j, then the simple loop a' := p'[xj, xj-] U /[x^, x^] must 
be non-contractible, and so separating some pair of trinodes of T from each other. In particular, 
by assumed connectivity of G', this implies that always > 1- Since at most 2k — 1 globally- 

shortest triedges of T span all the trinodes by Definition 9 a' (and consequently /[xj,Xj]) must 
cross at least one of them. Therefore, if a triedge p' repeats in S'/ at least 8/c^ times, then there is 
a globally-shortest triedge p of T such that p repeats in S/ at least i > 8k^/(2k — 1) > times. 

Let Y = (yi,... ,ye) (a subsequence of (x/,...,x™-^)) be the ordered sequence of points in 
which the arc of / intersects the arc of the triedge p. We say that an index i G {2, ...,£ — 1} is a 
switchback of Y if yi-i, yi+i both lie on the same side of yi on p. Up to symmetry, let the points 
on p be ordered such that yj+i lies between yi-i, yi. Since p is globally-shortest in T, we get (now 
regardless of contractibility of the induced loops) 


tw'{p[yi-i,yi]) < tw’{f[yi-i,yi]), 


and then 


tw'{f[yi-i,yi+i]) > 
> 


tw'{p[yi-i,yi]) +t^^'{f[yi,yi+i]) 
tw'{p[yi-i,yi+i]) +t^'{p[yi+i,yi]) +tw'{f[yi,yi+i]) 
tw'{p[yi-i,yi+i]) + 1 - 


Hence, if we locally re-route / along p[yj_i, yj+i], then we save the amount of at least I in the 
crossings of / with E(G). Note that this is not a contradiction to our choice of optimal drawing D 
yet since the change may introduce many new crossings of / with the rest of F. However, we cannot 
have more than ( 2 ) switchbacks in Y or we get a contradiction using Corollary (Observe that 
this usage of the corollary here is the only reason why we are restricting ourselves to unweighted 
sets F.) 

Since i > Ak^, there is a consecutive subsequence Y' CY of length F > i/(^^ — 1 > 8k without 
switchbacks. Without loss of generality, we assume Y' = (yi,..., y^/). Let y* := /[yi,yj+i] and 
g° := Qi U p[yi, yj_|_i], for i G {1,..., T — 1}. As argued before, each g° is a simple loop separating 
some pair of trinodes of T. Since no two edges of F cross more than once, there are at most k — 1 
indices i G {1,..., T — 1} such that gi is crossed by another edge(s) of F. 

Let X, y be the ends of the triedge p. Assume that we have i/yG{l,...,^' — 1} such that 
neither of gi,g°j separates x from y. Let Zi ^ ^ denote the set of trinodes of T that are separated 
by g° from x,y, and let Zj be defined analogously. We claim that Xi n Xj = 0. If not, then — 
up to symmetry — g° is separated from x,y by g°, except a possibly shared section of p[yi,yi+i]. 
The former is impossible by the Jordan curve theorem and the latter would mean that there is a 
switchback between i and j, which is again a contradiction. Since there are at most 2k —2 pairwise 
disjoint nonempty possibilities (e.g., singleton trinodes other than x,y) for the sets Xi,Xj, at most 
2k — 2 indices i G {1, — 1} are such that g° does not separate x from y. 
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Since > S/c, there exists a set of indices J C {1, — 2}, | J| > — 2(A; — l + 2fc —2) — 2 > 2/c, 

such that for every j £ J both the arcs gj,gj+i are not crossed by other edges of F and both 
gj,gj-^-l separate x from y. Let /o := f[yi-,yi'\ and po := p{yi-,ye]'i we get Y' C since there are no 
switchbacks in Y'. Observe also that g° H g°_^_l = {yj+i} since / is not self-intersecting and there 
is no switchback in Y'. Hence, up to symmetry, g° separates x from gj^i, and g°^-^ separates g° 
from y. It easily follows that gjUg°^^ forms the boundary of an arc-connected region of M^\(/oUpo) 
(a face of fo U po). Since at most 2fc — 2 of the faces of /o U po may contain a trinode of T other 
than X, y, there exists j & J such that, in addition to the above properties of J, the face cr bounded 
by g° U gj_^_i contains no trinode (see in Figure®. 

Our goal now is to re-route / along p[yj,yj^^ (i.e., “replacing” the part gj C /). Again, since p 
is globally-shortest in T, this move does not increase the number of crossings of / with E(G), and 
the T-sequence of / gets shorter. It remains to argue that we can avoid new crossings of / with 
F\ {/}. If any f'GF crosses p[yj,yj+i] then, since a contains no trinode, /' has to leave a as well, 
and the only possibility is across p[yj+i,yj+ 2 ] by the previous assumptions. Consequently, such /' 
can be re-routed along p[yj,yj^ 2 ], similarly to /, and no crossing with / is required (see again in 
Figure]^. Note, moreover, that even if two such edges f',f' G F cross each other in a, there is 
no problem and they will cross in their new routing in the same way. We have again reached a 
contradiction to our choice of D. □ 


3.3 Shortest routes in a sleeve, and crossing of routes 

Now, consider step |(IV)] of our outline—as mentioned before, we solve this step separately for each 
f = uv G F. To recapitulate, for trinodes n, u of a trinet T of G and a given proper T-sequence 
S from u to V, the task is to find a shortest route from u to u among those having the same 
T-sequence S. We cannot, in general, completely avoid repeating triedges in S and tricells in the 
sequence {Oq, 6i,..., 6m) in Definition 10 To prevent related technical difficulties, we use a similar 


workaround as in ; “lifting” the respective sequence of tricells into a universal cover as follows. 

Definition 13 (Sleeve of a T-sequence). Let {G',T) be a full trinet of a plane graph G, and 
consider a proper T-sequence S = {pi,p2, ■ ■ ■ ,Pm) from u to v determining the sequence of tricells 

For i = 0,1,..., m, let Li be a disjoint copy of the embedded subgraph 
Construct a plane graph L from the union Tq U • • • U Lm by identifying. 


( 6 * 0 , 01 ,..., 6jn) by Claim [Tl 
of G' U T induced by 0. 

for j = 1,..., m, the copy of the triedge pj in Tj_i with the copy of pj in Lj. We call L the sleeve 
of S in the trinet {G\ T), and we identify u and v with their copies in Lq and Lm, respectively. We 
make the unique face of L that is not covered by a copy of any tricell of T the outer face of L. o 


Observe that every route from tt to u in G' U T having its T-sequence equal to S can be easily 
lifted into a corresponding u-v route in the sleeve L of S. Conversely, any u-v route in L avoiding 
the outer face and crossing the copies of triedges in L at most once each, can be obviously projected 
down to G'uT to make a route with the T-sequence equal to S. In fact, we prove that some shortest 
u-v route in L must be of the latter kind, under the shortest-spanning property (cf. Dehnitionj^. 

Lemma 14. Let {G',T) be a shortest-spanning full trinet of an edge-weighted plane graph G, 
S a proper T-sequence between trinodes u,v of T, and let L be the sleeve of S. Let i be the length 
of a shortest route from u to v among those having the T-sequence S. Then, I is equal to the dual 
distance from u to v in L without the outer face. Furthermore, at least one of the routes of 
length i in L crosses the copy of each triedge from S in L exactly once. 
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Figure 3: Two consecutive crossings of the arc oi f = uv with a triedge p (which is locally-shortest) 
determine a contractible loop, and so / can be re-routed partly along p without inducing more 
crossings with G or with other edges of F. 


X 







Figure 4: The face bounded by g° U (C / Up) depicted in green: since there is no trinode 
in this face and neither gj,gj+i are crossed by other edges of F, it is possible to re-route / partly 
along p such that it avoids gj. Other possible F-edges entering the green face through a section of 
p must leave at the other end, and hence can be re-routed similarly to /. 
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Proof. Let tt be any shortest route from u to n in G' U T with the given T-sequence S. The copies 
of the faces and dual edges of vr lifted into the sleeve L give a route ttl from u to n which avoids 
the outer face of L. Obviously, the length of ttl equals the length of vr. 

Conversely, we aim to show that some shortest u-v route crosses the copy of each triedge of S 
in L exactly once. Assume a shortest route vri of length i from tt to n in L without the outer face. 
Recall from Definition 13 that L = Lq U ■ ■ ■ U Lm- For S = {pi,P 2 , ■ ■ ■ iPm), let {p'i,P 2 , ■ ■ ■ ^Pm) be 
the sequence of corresponding copies of the triedges of S in L, and let Pq := u, p^+i Note that 

each p', i G {1,..., m}, connects two vertices of the outer face of L, and so p' separates p[_i from 
p'i^i- In particular, every u-v route in L which avoids the outer face must cross each of p'^,... ,p^. 

Let i be the maximum index such that p^ is crossed by vri more than once. Then there is a 
subpath fji C TTi stretching between two consecutive crossings of vri with p' and contained in Lj. 
We turn vri into 7r2 by re-routing the subpath ai along the boundary p' in Lj_i. Since p* is locally- 
shortest in the trinet T, the length of 7r2 equals the length of tti. By induction on the number of 
excess crossings of 7r2 with copies of the triedges, we can then get a u-v route ttq of length £ such 
that TTo crosses the copy of each triedge from S' in L exactly once. 

Finally, the route ttq projects down to a route of length i from u to v in G' L) T having the 
T-sequence S. □ 


Regarding the shortest path computation, we note that if the edge-weights are given in unary or 
are bounded by a constant, a simple adaption of BFS achieves the job. Otherwise, we can use the 
algorithm of Klein et al. 27 since L is planar, or Thorup’s algorithms |31| since we have integral 


weights. Altogether we obtain: 


Corollary 15. Let {G',T) be a shortest-spanning full trinet of an edge-weighted plane graph G, 
and S a proper T-sequence between trinodes u,v ofT. A shortest u-v route among those having the 
T-sequence S can be found in the geometric dual graph of the sleeve L of S inO{\S\ • | A'(T)| • |1/(G)|) 
time, using a linear time shortest path algorithm. 


Observe that in our case the term |S| • |iV(T)| is bounded by a function ot k = \F\. 


Proof. By Definition the size of the full trinet {G',T) is C>(|A^(T)| • |K(G)|), and the sleeve is 
composed of |S| copies of subgraphs of G' U T, and so C>(|5| • |A^(T)| • |1/(G)|) bounds the size of 
the sleeve L. Therefore, any linear time shortest path algorithm applicable in our situation, e.g. 
the aforementioned Klein et al. pTl or Thorup 31 , does the job. □ 


Finally, it remains to address step |(V)[ Consider a 4-tuple of distinct trinodes u,v,u',v'. Let tt 
be a u-v route and tt' be a u'-v' route. We say that an arc b follows the route vr if 6 is contained 
in the alley of tt and b intersects the faces forming the alley exactly in the order given by vr (recall 
that a route is technically a dual walk and hence, possibly, some face might repeat in vr). We say 
that the pair of routes vr, vr' is non-crossing, if there exist a u-v arc b following vr and a u'-v' arc b' 
following vr' such that 5 n 6' = 0. In order to characterize possible non-crossing pairs of routes in 
terms of their T-sequences, we bring the following definition: 


Definition 16 (Crossing certificate). Let {G',T) be a full trinet of a plane graph G, and let vr be 
a route from n to n and vr' be a route from u' to v' in G' U T, where u, v, u', v' are distinct trinodes 


,p^) of vr' are proper and 
For technical reasons, let 


of T. Assume the T-sequences S = {pi,... ,pn) of vr and S' = {p'^,. 
let {Oq, ..., 9n) and (0g ,... ,6'^) be their tricell sequences by Claim 11 
po := u, pn+i ■■= V and p'g := u', := v'. 

A crossing certificate for S, S' is a triple of indices (c, d, m) where c,d,m > 0, c -|- m < n, 
d -\- m < I, such that the following holds: 
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a) Oc+j = for 0 < j < m, but pc + p'^ and Pc+m+i + p'd+m+n 

b) the triple pc,p'^-,Pc+i occurs around the tricell 9c = 9'^ in the same cyclic orientation as the 
triple Pc+m+i,p'd+m+i^Pc+m occurs around 9c+m = G'd+m- 

Furthermore, a crossing certificate for the same sequence S and the reversal of S' from v' to u' is 
also called a crossing certificate for S, S'. o 

Dehnition[^ deserves a closer explanation. Assume that a crossing certificate satisfies 0 < c < n 
and 0 < d < Then all four elements pc-,p'diVc+i-,p'd+i triedges of the same tricell 9c = 9'^, and 
since Pc+i 7 ^ Pc p'd / P'd+i-i Pc+i = p'd+v Hence m > 0 and the situation is such that S 

and S' “merge” at 9c where (up to symmetry) S comes on the left of S', and they again “split” at 
9c+m where S leaves on the right of S', thereby “crossing it”. The full definition, though, covers 
also the boundary cases of crossing certificates for which c G { 0 ,n} or d G { 0 ,^} (or both), and 
when S and S' may have no triedge in common; those can be easily examined case by case. 


Lemma 17. Let {G', T) he a full trinet of an edge-weighted plane graph G, and Ui, Vi, i = 1,2, be 
four distinct trinodes. Assume that Si from Ui to Vi are proper T-sequences. In G' U T, for i = 1,2, 
there exist routes iTi from Ui to Vi having the the T-sequence Si, such that vri,7r2 are non-crossing, 
if and only if there exists no crossing certificate for 5i, ^ 2 . 


Suppose we have two proper T-sequences S,S' as in Definition 16 Referring to Definition 16 


point a), we call the tricells 9c+j = ^'d+j 0 < j < m the central tricells of the crossing certificate 

(c, d, m). 


Proof. Let Li, i = 1,2, be the sleeves of Si in (G', T). Assume that (c, d, m) is a crossing certificate 
for Si, S 2 , and let R = {9c, ■.., 9c+m) be the sequence of the central tricells of this certihcate. Let 
Ki 'L Li, i = 1,2, he the plane subgraphs consisting of the copies of the tricells from R in the sleeve 
Li. Note that R may repeat the same tricell several times, but in Li we have got independent 
copies of the possibly repeated tricells. We may also assume that Ki = K 2 since they are both 
made of copies of the same sequence R of tricells. 

says that (c, d, m) is a crossing certificate iff the elements 


In the above view. Definition 16 


Pc,p'djPc+m+i,Pd+m+i appear on the outer face of Ki = K 2 in this cyclic order. Hence, by Jordan’s 
curve theorem, if there is a crossing certihcate for Si,S 2 , then ' 7 ri, 7 r 2 cannot be non-crossing. 


Conversely, we show how to build non-crossing vri, 7 r 2 if there is no crossing certihcate for Si,S 2 . 
For each tricell 9 of T, bounded by triedges qi,q 2 ,Q 3 , we choose arbitrary three edges Cj from qj, 
j = 1, 2, 3, and arbitrary three internally disjoint dual paths cri, 2 , < 71 , 3 , ( 72,3 contained in 9 such that 
iTjj is a dual path in G' U T connecting the face incident with e* to the face incident with ej. 
Furthermore, we denote by xj the trinode of 9 opposite to qj and we choose another three arbitrary 
dual paths pi, P 2 ,P 3 contained in 9 such that pj is a dual path in G' U T connecting a face incident 
with Xj to the face incident with Cj. We call chosen cri^ 2 , < 71 ^ 3 , 1 x 2 , 3 , Pi, P 2 , P 3 the representative dual 
paths of the tricell 9. 

For the proper T-sequence Si, i = 1,2, we simply compose the route vTj of the apropriate 
representative dual paths of the tricells determined by Si. It is routine to verify that these vri, 7 r 2 
are non-crossing, if and only if there exists no crossing certificate for S'!, ^ 2 . □ 


3.4 Summary of the algorithm 

We are now ready to put all of the above results together, in order to summarize the overall algo¬ 
rithm to solve r-MEI, see AlgorithmBased thereon, together with Lemmas 12,^ Corollary |15[ 
and Lemma ITtI we obtain: 
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Theorem 18. Let G be a conneeted plane graph with edge weights w: E{G) —)• N+ U {oo}, and F 
a set of new edges (vertex pairs, in fact) such that w{f) = 1 for all f ^ F. Algorithm^ finds an 
optimal solution to w-weightedT-MEl{G,F), if a finite solution exists, in time . 


Input: 

Output: 

( 1 ) 


( 2 ) 


(3) 


(4) 

(5) 


a plane graph G, edge weights w: E[G) —>■ N+ U {oo}, new edge set F s.t. w{f) = 1 for / G F. 
an optimal solution to (w-weighted) r-ME^G, F). 

Compute a full trinet {G',T), N{T) := V{F), with the shortest-spanning property of T. 

(a) Pick any trinode x G N{T) and greedily compute globally-shortest triedges (Def. from 
X to all other trinodes, using a simple shortest path computation. 

(b) The remaining triedges can be greedily computed as locally-shortest, one after another. 
For each f = uv G F: 


(a) Compute §/ as the set of all of its possible proper T-sequences from u to v. The size of §/ 


12 


b) as 


is bounded due to Lemma 
(b) For each S' G §/, compute a shortest u-v route ns in G' LI T among those having the 


F-sequence S (where the length function is induced by w), using Corollary 15 
For each possible set IP = {S/}/g_F with S/ G §/: 


(a) Check, for each pair /,/' G F, whether there exists a crossing certificate for Sf,Sfi 
(e.g., using brute force by Def. [T^. 

Let Xy be the set of pairs {/,/'} for which such a certificate has been found. 

(b) If any pair {/,/'} G X‘j> requires more than a single crossing (which can be found by 
checking again for two “independent” crossing certificates of S/,S/'), let cry := oo. 

(c) Otherwise, let cry := |Xy| -P where len^ denotes the length function in 

the geometric dual of G' U T induced by w. 


Among all IP considered in (3) 


pick the one with smallest cry < oo. Let this be 1P° = 


In the plane graph G, realize each edge / G F following its respective route ns °, such that the 
overall resulting weighted number of crossings is cryo. 


(a) (By minimality, no ns° will be self-intersecting.) Using well-known postprocessing— 
removing consecutive crossings between /, f' by re-routing f' partially along / or vice 
versa—allows to avoid multiple crossings in pairs from Ay and to make remaining pairs 
from F crossing-free. 


Algorithm 1: Algorithm to solve the (weighted) rigid MEI problem. 


Before giving the proof, we need a deeper understanding of the concept of non-crossing routes and 
crossing certificates, and a detailed specification of the step (3( (b) of Algorithm 


By adapting the arguments of Lemma pG} one can actually get the following slight strengthening: 


Claim 19. Let {G',T) be a full trinet of an edge-weighted plane graph G, and Ui,Vi, i = 1,2, be 
four distinct trinodes. Let ni, i = 1,2, be a ut^i route in G' LiT. If there exist simple Ui^i arcs bi, 
i = 1,2, following the route ni, such that bi intersects 62 in exactly one point x and they properly 
eross in x, then there exists a crossing certificate for the T-sequences of ni and n 2 . □ 


We say that there exist two independent crossing certificates for the T-sequences S,S' if there 
are crossing certificates {c,d,m) and {c',d',m') for S,S' (each one up to possible reversal of S' as 
in Definition 16), such that the set of central tricells of {c,d,m) is disjoint from the set of central 


tricells of {c', d', m'). The following can then be straightforwardly obtained from Lemma 17 


Claim 20. Let {G',T) be a full trinet of an edge-weighted plane graph G, and Ui,Vi, i = 1,2, be 
four distinct trinodes. Assume that Si from ui to Vi are proper T-sequenees. In G' U T, there exist 
simple arcs bi from Ui to Vi such that, for i = 1,2, 
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- bi is contained in the alley of a Ui^i route having the T-sequence Si, and 

- bi intersects 62 in at most one point, 


if and only if there exists no two independent crossing certificates for Si, S 2 - 


□ 


The implementation of step (3^ b), using Claim 20, hence simply checks by brute force for the 
existence of two independent crossing certificates for Sf,Sf>. 


Proof of Theorem 18. In this proof, we will use some of the terminology and notation from the 


proof of Lemma and refer to the notation of Algorithm 

Consider arbitrary T = as in the step|(3)[ The value of cry computed in the step ^ 'c) 


provided that cryo < 00 , is a lower bound on the number of crossings of any feasible solution of 
r-MEI(G, F) such that, for each f £ F, the T-sequence of the arc / is exactly Sf. This fact follows 
directly from Lemma 17 (for the part |Ay|) and from Lemma 
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(for the part 

By Lemma 12, there is an optimal feasible solution D to r-MEI(G, T) such that, for every 


f £ F, the arc of / in T has its T-sequence (with respect to the trinet T from step (1)) equal 
to some proper Sf £ §/ as computed in ( 2 ), and the step ^'b) does not apply to these values 


by Claim 20 Consequently, cryo < r-ins(G, F) by the lower-bound argument from the previous 
paragraph. Hence if we can prove that the step |(5)| indeed can compute a drawing D° of G + F 
with cryo weighted crossings, provided cryo < 00 , then we complete the proof of Theorem 18 


For / £ F, let bf denote a realization of the edge / as an arc following the route of tts° 

By the minimality 


(such that the T-sequence of bf is Sj), before the postprocessing step (5 


choice in step (4), we can be sure that bf does not cross itself: the self-crossing would induce 


a non-contractible loop with at least one crossing over G', but then there exists a T-sequence 
S'j, C Sj —the sequence SJ without the triedges forcing the loop. Replacing these two sequences in 
1 P° results in a smaller number of crossings on / while not increasing the number of crossings at 
any term in the summation considered in step (3)|[c) 


Let Di denote the drawing of G-l-T made of G and {bf : f £ F}. Observe that cr^i^ {E{G),F) = 

'Ylf^F 

Fix some f,f' £ F, f = uv and f = u'v', such that bf,bft cross each other (properly) more 
than once. Consider the point set p := bfU bf/ with the outer face incident with, say, the trinode u. 
If any of the bounded faces of p contained a trinode (which cannot be u) of T, then we could “split” 
the T-sequences SJ and SJ, into SpSj and Sj,,S‘j, each, such that for each of the pairs Sj,Sj, 
and S‘f,Sj, there would exist a crossing certificate by Claim pT This would, in turn, provide two 
independent crossing certificates for the T-sequences S^SJ, £ and by the step (3^ b), it would 
contradict cryo < 00. 

Consequently, all the bounded faces of p = bfU bf/ are free of trinodes of T. This, in particular, 

means that if we construct from 6/ by re-routing it along a section of bf/ between two consecutive 

shared points of bfCibf/, then the T-sequence of would again be SJ. Moreover, since bf, bf/ have 

been chosen from their respective shortest routes, the t-weight of b\ would be equal to lenw{'Fs°) 

•1 f 

(which is the t-weight of original bf). Iterating this process, we arrive at a drawing D 2 of G + F 
satisfying the following: 

i. no two edges of F in D 2 cross more than once, 

ii. crD,{E{G),F)<ciDAE{G),F). 

What remains is to observe that two edges /, /' £ F properly cross each other in D 2 only if 
ifjf'} £ Indeed, if /,/' properly cross in D 2 , then this crossing is the only one and there 

exists a crossing certificate for vr^o , fs°^ by Claim 
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and then {/, /'} £ Ayo due to the step (3 
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To summarize, we get 


cryo < j:-ms{G,F) < ci{D 2 ) = ctd 2 {E{G),F) + ci:D 2 iF,F) 

< CVDi {E{G),F) + \Xyo I = cryo 

which proves optimality of the solution computed by Algorithm 

Finally, we discuss the runtime bound of Algorithmic Let k = |T|. Step (1) is performed 
in time 0{k‘^ ■ |y((S')|) using 3{2k) — 6 calls to a lin ear shortest path algorithm. Step (2) takes 
time 0{k ■ • k\V{G)\) by Corollary 15, Step (3) is iterated times, and each 


iteration takes time polynomial in k (independently of G) even by brute force. Step |(4)| takes only 
time 0 { 2 P°^y^^k^y Finally, step (5) performs k computations in 0(\V{G)\), to realize each f ^ F 


in G, and then a number of concurrent re-routings which can be bounded by an amortized analysis: 
every of the k routes is of length Oi\V{G)\) and each element of it could be re-routed at most once 
towards each of the k — 1 remaining routes, summing to 0{k‘^ ■ |C(G)|). 

The above analysis sums up to overall 0(2"^°^^^^^ ■ |C(G)|) time. □ 


4 General MEI 


Now, we may turn our attention to the general MEI problem, where the embedding of the planar 
graph G is not prespecified. See also the appendix for details. Recall that triconnected planar graphs 
have a unique embedding (up to mirroring), but already biconnected graphs have an exponential 
number of embeddings in general. As it is commonly done in insertion problem since [^, we will 
use the SPR-tree datastructure (sometimes also known as SPQR-tree) to encode and work with 
all these possible embeddings, 
work of 


It was first defined in slightly different form in 
It can be constructed in linear time 


15 , based on prior 


19,25 and only requires linear space. 


Definition 21 (SPR-tree, cf. jC). Let G be a biconnected graph with at least three vertices. The 
SPR-tree T of G is the unique smallest tree satisfying the following properties: 


i) Each node v in 7 holds a specific (small) graph = {Vy^Ey), with Vy C V{G), called a 
skeleton. Each edge e of Ey is either a real edge e G E{G), or a virtual edge e = xy ^ E{G) 
(while still, x,?/ G V{G)). 

ii) T has three different node types with the following skeleton structures: (S) Sy is a simple 
cycle; (P) Sy consists of two vertices and at least three multiple edges between them; (R) 
Sy is a simple triconnected graph on at least four vertices. 


iii) Eor every edge uy in T we have 114 n V)j| = 2. These two common vertices, say x,y, form a 
vertex 2-cut (a split pair) in G. Skeleton Sy contains a specific virtual edge G E{Sy) that 
represents the node p and, symmetrically, some specific Cy G E(5^) represents ly; both 
have the ends x, y. 


iv) The original graph G can be obtained by recursively applying the following operation of 
merging: Eor an edge up G E{7), let e^, Cy be the pair of virtual edges as in (iii) connecting 


the same x,y. A merged graph (Sy U 5^) — {e^,ey} is obtained by gluing the two skeletons 
together at x,y and removing efj,,ey. o 

The central theorem of states that we can find an optimal embedding to insert a single edge 
uv by looking at the shortest path in 7 between a node whose skeleton contains u and a node whose 


17 












skeleton contains v. For each skeleton along this path, one considers the partial routes between the 
virtual edge representing u (or u itself) and the virtual edge representing v (or v itself). In case 
of S- and P-nodes this route requires no crossings (by choosing a suitable embedding in the latter 
case); for an R-node the route is a shortest path in the dual of its skeleton: if the primal edge 
is an original edge, the length of its dual edge is the primal edge’s weight; if the primal edge is 
a virtual edge xy, representing node y, the length of its dual edge is the minimum-xy-cut in P^, 
where we is the pertinent graph of y arising from merging all skeletons of the subtree rooted 
at //, minus the edge e^,- By picking any embedding of and computing a shortest dual path 
through it, we can compute this cut size in linear time. See for details. 

We consider our SPR-tree T of G rooted at any node, and devise a dynamic programming 
scheme to solve MEI bottom-up over 7. We observe that every non-root skeleton Si, contains a 
virtual edge that represents its father in 7. Any further virtual edges correspond to children of i' 
in T. Since we already know how to solve r-MEI, it shall suffice to describe which r-MEI problems 
we need to solve at each SPR-tree node n (including the root node), assuming we already solved 
the corresponding subproblems at their children. The overall MEI solution can then be obtained 
by selecting a solution in the root with the least number of crossings. 

We say a virtual edge in Si, is dirty if it contains an end vertex of a new edge f G F. Hence, at 
most 2k edges (a constant number) of Si, are dirty. For R- and S-nodes, we only have to consider 
their unique (up to mirroring, in case of R) embeddings. A P-node whose skeleton contains p edges, 
however, allows {p — 1)! embeddings. Based on the following claim (which can be shown with a 
straight-forward redrawing argument), we only need to consider up to (2A:)! embeddings for each 
P-node, which is constant for constant k. 

Claim 22. Let v be a P-node in the SPR-tree 7 of G. There is an optimal embedding of G for the 
MEI problem, where all non-dirty virtual edges are consecutive in the embedding of S^. 

Let S^ be a considered embedding of Sy. Consider each virtual edge xy, representing node 
pL (possibly, fi is the father node), in Sy. If xy is not dirty, we set its weight to the size of the 
minimum-xy-cut in P^ (as for the single edge insertion case). If xy = e is dirty, we modify it with 
the following gadget: Set the weight of e to oo, and add two new side edges e',e" connecting x 
and y. One is directly to the left, the other directly to the right of e. We will further modify these 
side edges in the following. 

Consider what can happen at the subdrawing of (embedded) component P^ in the context of 
whole G when considering any specific new edge f G F: (i) if / has exactly one end in P^, it will 
enter the component; (ii) if / has no end in P^, it may cross through the component; (hi) if / has 
both ends in P^, it may leave the component and re-enter P^ at another position. Furthermore, and 
in contrast to the single edge insertion, it may happen that / (independent of its end points) crosses 


that the latter number is bounded by a constant, depending only on k. Hence, there are only a 
bounded number of enterings/leavings at P^, and we can simply consider all such possible situations 
(including all possible orders of the enterings/leavings). For each such situation, we now subdivide 
the edges e' and e" accordingly: chiefly put, if, e.g., we consider the case of an edge f = uv coming 
from a vertex u 0 V{P^) and crossing P^ twice before finally entering it to reach v G V{P^), we 
generate (for /) overall five vertices on e' ,e" , say vi,. .. ,v^. Within the context of the dynamic 
programming subproblems at P^, we then consider (for each embedding of S^) the r-MEI problem 
w.r.t. subedges viV2, U3U4, v^v. Overall, we have to store the best solution (over all embeddings of 
S^) for each r-MEI problem constructed of all such subedges, for each edge f G F, each possible 
number of crossings of / through P^, each possible assignment of thereby induced subdivision 
vertices to e' or e" , and all possible orders at e' and e" . Within a subproblem at v, we then 


P^ multiple times. However, since we consider a fixed embedding S'(,, we know from Lemma 
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consider the edge uvi instead of / (in fact, this edge may be further split into several subedges 
due to further dirty virtual edges in Su considered to be crossed by /, and/or if u is contained in 
a pertinent graph of another virtual edge). 

Hence, we only need to store a constant (bounded by a function in k) number of solutions 
at each SPR-tree node. Each solution can be obtained using the above algorithm for r-MEI in 
0{\V{G)\) time, and there are at most 0{\V{G)\) SPR-tree nodes. Instead of the naive quadratic 
runtime bound, we even achieve a linear runtime bound by observing that the union of all skeletons 
is still only of linear size. We obtain, as given in the introduction: 


Theorem 23 (The biconnected case of Theorem [^. Let G be a planar bieonnected graph on n 
vertiees, and F a set of k new edges (vertex pairs, in fact) where k is a constant. We can solve 
MEI(G, F) in 0{n) time. 


Por essentially all known insertion algorithms (in particular single edge insertion 20 
insertion [^, and MEI approximation 


vertex 

[10] ), one can typically first describe the case of biconnected 
graphs (using SPR-trees). Then, it is relatively straight-forward to lift the algorithms to connected 
graphs, by considering BC-trees (see below). Interestingly, this seems much more complicated in 
case of exact MEI: 

Consider the well-known block-cut tree (BC-tree) to decompose any connected graph into its 
blocks (biconnected components). Using analogous techniques as in 10 , we extend our dynamic 


programming approach by amalgating the BC-tree with the blocks’ respective SPR-trees, to obtain 
a linear-sized con-tree, with an additional node type C, for cut vertices. In our bottom-up approach, 
at a cut vertex c, we need to consider all possibilities to “glue” the c-incident dirty blocks (blocks 
with at least one end of some edge f € F) together. However, we cannot easily bound this number 
by a function purely in k: we not only have to consider all orders of these blocks, but also all 
possible nestings, which introduces a dependency on Acutj the maximum degree of the cut vertices 
in G. Hence, for only connected G, we obtain the slightly weaker result: 


Theorem 24 (The connected case of Theorem [^. Let G be a planar connected graph on n vertices, 
and F a set of k new edges (vertex pairs, in fact), where k and the maximum degree of the cut 
vertices in G are constant. We can solve MEI(G, F) in 0{n) time. 


As sketched in the main body of the paper, we hrst develop a dynamic programing algorithm 
over the SPR-tree decomposition T of planar G. This algorithm considers dirty nodes bottom-up; 
a decomposition node v is dirty if its pertinent graph contains at least one vertex incident to F. 
Observe that if a node is dirty, so is its parent. The root node (whose pertinent graph we may 
define as G itself) is always dirty. 


Subproblems at non-root nodes. We start with formally defining the subproblems to be solved 
and stored at each dirty non-root decomposition node. Let v be such a node and let e = xy G E{Su) 
be the virtual edge in the skeleton of v corresponding to its parent node q. Recall that the pertinent 
graph Py arises from Sy by merging the skeletons of the subtree rooted at v and removing the sole 
remaining virtual edge (e). We consider the 3-partition of F into Fq, Fi, F 2 , where Fq are the edges 
without an end in V{Py) \ {x,y}, Fi are the edges with one end in V{Py) \ {x,y} and the other 
not in V{Py), and F 2 are the edges with one end in V(Py) \ {x,y} and the other in V{Py). 

By definition, the graph F+ -.= Py + e is planar, and e represents the “rest of the graph” disjoint 
from Py. We are, intuitively, interested in the best embedding Pf of F+ to 

(a) route the edges of Fi from a side of e to its end in V{Py) \ {x,y}; observe that we may care 
from which side of e the new edge emanates. 
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But these are not the only routes to consider in an optimal solution: 

(b) edges uv G F 2 may be routed completely within Py, or go from u to some side of e (into the 
“rest of the graph”), and from some side e (from the “rest of the graph”; either the same or 
the other side) to v\ 

(c) any edge of F may be routed through i.e., from one side of e to the other side, without 
crossing e. 


Formally, we can define a routing query as a pair (s,t), where s and t are each either referencing a 
specific side of e or a vertex in V{Py) \ {x,y}. We will use e',e" to denote the two different sides 
of e. In such a routing query, we ask for a routing of a new edge between s and t in P+, without 
crossing over e. 

Lemma 12 (which holds for every hxed embedding, and hence for each possible embedding) 
showed that a triedge of trinet T is crossed at most 8 /c^ = 0{poly{k)) times. When computing 
a shortest route (w.r.t. some T-sequence) between two succeeding triedges, we clearly have the 
property that any edge within the corresponding tricell is crossed at most once. Hence: 


Corollary 25. In an optimal solution to MEI(G, P), each edge f £ F crosses any edge e £ G at 
most ^ := poly{k) times. 


Since this corollary also holds for virtual edges in a skeleton, we have the same upper bound 
for crossings through a two-connected component P^. 

In our dynamic programming scheme, we will hence—for each possible set of routing queries— 
store the minimum number of crossings necessary over all embeddings of P+. A specihc set of 
routing queries (to be described in details below) is hence a subproblem, and the corresponding 
number of crossings (together with the embedding of P+ and the corresponding routings, if desired) 
is a subsolution. It remains to discuss the number of subproblems for u. 


Lemma 26. Each subproblem specifies at most O{poly{k)) routing queries. The total number of 
subproblems to consider at any node v is bounded by O{poly{k)\). 

Proof. Consider the routing types (a)-(c) as above: 

(a) For each edge f = uv £ Fi with u £ V[Py) \ {x,y}, we have to pick one of the two routing 

queries {e”,u). 

(b) For each edge f = uv £ F 2 , we have to pick one out of five options: (i) a single routing query 

{u,v); (ii)-(v) two routing queries (n,e^^^), {e^‘^\v), with £ {e',e"}. 

(c) Finally, for each / G P—except for those P 2 -edges that picked option (i)—we have additional 
up to ^ routing queries. Each such additional query is of one of four types: (e^-^^e^^^), with 

g(l), e(2) g |g/^ gWj^ 

Overall, this gives up to r := |Pi| -|- 2|P2| + ^\F\ = O{poly{k)) routing queries. 

The number of choices for such a set of routing queries is at most • (4^)l^l = 

However, up to now we did not consider a crucial interplay of these individual routing queries: We 
need to take all possible orderings of the edges emanating from a side of e into account: Sides of 
e arise at most 2r times over all queries, and we hence have at most (2r)! orderings to consider. 
Thus, we overall obtain C)( 5 P°^y(C . poly(k)\) = O{poly{k)\) subproblems. □ 
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Dynamic programming and root node. Finally, we have to describe how to use these sub¬ 
problems to efficiently compute MEL The validity of this approach for non-dirty pertinent graphs 
was already established in 20 . As mentioned, we consider dirty nodes bottom-up. 

Let V be the considered SPR-tree node with skeleton Sv Let G E{Su) be the virtual edge 
corresponding to i/’s father g (if it exists), and ei,..., (e'^,..., e^,) the dirty (non-dirty) virtual 

edges in corresponding to the children ^i,... respectively). We need to show 

that we can solve each subproblem at v purely using Si, and the solutions to the subproblems of 
the dirty children. In particular, we may not expand the skeleton to the pertinent graph (for which 
the i^-subproblems are actually defined). 


Subproblems, embeddings, and the root. Assume, v \s a, non-root node, then we have to solve 
X := O{poly{k)\) many subproblems. For each subproblem, we are given a set of r := O{poly{k)) 
routing queries, and want to find the optimal solution over all embeddings of P+. As a first step, 
we recall that there are only a bounded number of embeddings for Sy (1, 2, and (2A;)! in case of an 
S-, R-, and P-node; let x' = C1((2A:)!)), and we may hence enumerate each one explicitly. 

The routing queries of the considered subproblem give rise to the following gadget: Set the 
weight of eg to oo, and introduce two edges e' and e" parallel to eg, one directly to its left, one 
directly to its right. Now subdivide these to edges such that there is a vertex on e' (e") if a routing 
query specifies the edge side e' {e”). Furthermore, these vertices are ordered according to the 
specification of the subproblem. Let S'l, denote the embedded graph arising from this construction; 
we do not consider eg as a virtual edge in the following any more. Instead of considering the original 
new edge set F, we will now consider the routing queries (s, t) as edges st G F' (a new set F') to 
be inserted. 

If n is the root node, we also have to consider all its skeleton’s possible embeddings indi¬ 
vidually, but there is no specific subproblem to consider and we simply set F' := F without any 
gadget construction. From now on, we do the same steps, independent on whether considering the 
root node, or a specific subproblem at a non-root node. 

Virtual edges. For each non-dirty virtual edge e' = a6, 1 < z < /, we set the weight of e' to the 
minimum-aft-cut in the pertinent graph of //(. Note that these values can be constructed bottom 
up in overall linear time as a preprocessing. 

Now, for each dirty virtual edge ei, 1 < i < i {(. < 2k), we construct a gadget analogous to the 
gadget for eg: Edge Cj gets weight oo, we add two edges e(,e" left and right of Cj, and subdivide 
them according to a subproblem at /Zj.—To do this, we have to enumerate all possible choices of 
subproblems at all virtual edges. So this construction yields x" ■= O{{poly{k)\)‘^^) different choices, 
each of which we consider individually. Observe: If, for some edge f € F that resides within some 
Pfij, we chose a routing query of type (b) (i), we do not consider subsolutions at any virtual edge 
where there are type (c) queries w.r.t. f. We call such an edge f a suppressed edge. 

We denote the so-modified plane graph by S”, and now have to decide what happens to our 
new edges F'. Each edge in F' corresponds to some edge in F. Furthermore, we add each non- 
suppressed edge / G F to E' if it has no corresponding edge in F'. We observe that for each vertex 
w G Py \ Sy that is an end in F', there is a unique replacement vertex r{w) in Sy —it arises from a 
query {r,r{w)) (unoriented) within a subproblem at some dirty virtual edge. 

For each original edge f = uv G F, we hence get a partial order of routing queries corresponding 
to it: either u (v) or its replacement vertex r{u) {r{v), respectively) is in S^, so we start (end) there. 
There may or may not be a routing query starting at u (ending at v), which we would update to use 
r{u) {r{v)) instead of u (v). Now, between this start and end, / may have to “visit” former queries 
of type (c) (whose ends are now represented by subdivision vertices at edges e'^e'-, 1 < j < tj. 
While these former queries are totally ordered for each individual dirty virtual edge, it is unclear in 
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which order / visits the different virtual edges. We will enumerate all possible orders to visit each 
of the 0{k) dirty virtual edges up to 0{^) times; there are hence x”' ■= O = O{poly{k)\) 

different visit orderings for each edge of f € F. Every visit order induces an unambiguous set Tj 
of (new) routing queries to draw part of / within Si^: from /’s start to the vertex representing the 
beginning of a former query, from the vertex representing the end of the last former query to the 
beginning of the next former query, and so on, until finally from the vertex representing the end of 
the last former query to /’s end. Such a set Tf hence has size at most 0{k^) = O{poly{k)). 

So, finally, we obtain an instance r-ME^S'^,, F"), where |E(S'(,)| = 0{\Si,\-k^) = 0{\Su \ ■poly{k)) 
and F" is the set of all routing queries (interpreted as unordered new edges) obtained from F' by 
considering each Tf (for all / G F). We have |F"| = 0{kpoly{k)) = O{poly{k)). The total cost 
of the considered subsolution (and also for the solution at the root node) is the minimum number 
of crossings over all possible r-MEI instances constructed as above plus the numbers of crossings 
given by the corresponding individual subsolutions realized at the dirty virtual edges. We have: 

Lemma 27. We settle the root node—and any speeific subproblem at a non-root node—with 0{x' ■ 
x" • x'"^) calls to r-MEI. We settle each dirty non-root node with 0(x ■ x'' x" ' x"'^) calls to r-MEI. 

Theorem 28 (Detailed version of Theorem |23[). Let G be a planar biconneeted graph on n vertices, 
and F a set of k new edges (vertex pairs, in fact) where k is a constant. We ean solve MEI(G, F) 
in 0{n ■ {poly{k)\)^^^'>) = 0{n) • time. 

Proof. Eirst, due to Theorem |18t each individual r-MEI instance in our setting can be computed 
within 0{\V{Si,)\ ■ poly{k) ■ 2^°^^)) = 0{\V{Si,)\ ■ time. Eurthermore, the union over all 

SPR-tree skeletons has still linear size 0{n). We hence obtain the overall runtime 

0{x ■ x! ■ x!' ■ x'"" • n • = 0{n • poly{k)\ • {2k)\ • (poly{kf)^'^ • (po/|/(A:)!)'= • 

= 0{n ■ {poly{k)\)^'^^^) 


□ 

Connected Case. Until now, we only considered biconneeted G. In case of only connected G, 
we can first decompose (in linear time) G into its biconneeted components (blocks), and establish a 
BC-tree 23. This tree has two types of nodes: Eor each block of G, we have a node of type (B); for 
each cut vertex in G, we have a node of type (C). We have an edge /Sy in 23 if, and only if, /3 is a 
B-node, 7 is a C-node, and the block of /3 contains the cut vertex of 7 . We may root 23 arbitrarily 
at any dirty block; we say a block is dirty if it contains at least one end of F (other than possibly 
its parent cut vertex). Clearly, we can iteratively prune non-dirty B-leaves. 

Now, we can construct a combined tree C: Eor each block B in G, we construct (and root) its 
SPR-tree T^. In 23, we replace each B-node with the root vertex of the block’s corresponding SPR- 
tree. Now, we can run the dynamic programming algorithm over S instead of a single SPR-tree 
tree. 

Let be a non-C-node whose parent is a C-node 7 corresponding to cut vertex c G S^,. We need 
to redefine the subproblems to consider at n: instead of considering routing queries that attach 
to one of the two sides of the parent virtual edge, our routing queries may now attach to c in a 
specified order and through specified faces incident to c. We therefore introduce the gadget—for 
each considered embedding of —obtained by planarly replacing c by a simple cycle G. The 
c-incident edges are attached to G such that the contraction of G again gives S^,. When considering 
the routing queries, instead of the two choices of the side of the parent virtual edge, we now hence 
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have a 5(c) = C>(Acut)-fold choice over the segment of C where to attach to, where Acut denotes 
the maximum degree over all cut vertices. 

In our dynamic programming, we will perform no operation at C-nodes, but let u now be a 
node with a C-child 7 corresponding to cut vertex c G Si,. Analogous to above—in each considered 
embedding of —we planarly replace c by a cycles C. On C, we realize all subsolutions of 
all (at most 2k) children of 7, in all possible combinations. Except for these modifications, the 
algorithm remains unchanged, and we obtain: 

Theorem 29 (Detailed version of Theorem |24[) . Let G he a planar connected graph on n vertices, 
and F a set of k new edges (vertex pairs, in fact), where k and Acut— the maximum degree of the 
cut vertices in G—are constant. We can solve MEI(G, E) in 0{n ■ Acut®^^^ • {poly{k)\)^^^'l) time. 
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